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ABSTRACT 


A prediction of the three-dimensional turbulent flow in the NASA 
Low-Speed Centrifugal Compressor Impeller has been made. The calcula- 
tion was made for the compressor design conditions with the specified 
uniform tip clearance gap. The predicted performance is significantly 
worse than that predicted in the NASA design study. This is explained 
by the high tip leakage flow in the present calculation and by the 
different model adopted for tip leakage flow mixing. The calculation 
gives an accumulation of high losses in the shroud/pressure-side quad- 
rant near the exit of the impeller. It also predicts a region of merid- 
ional backflow near the shroud wall. Both of these flow features should 
be extensive enough in the NASA impeller to allow detailed flow measure- 
ments, leading to improved flow modelling. Recommendations are made for 
future flow studies in the NASA impeller. 
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NOMENCLATURE 


A,B,C 

grid indices in throughflow, blade-to blade, and hub-to-shroud 
directions 


C P 

specific heat capacity at constant pressure 


D 

diameter 


h 

height of blade tip clearance gap 


H* 

rothalpy 


AH 

change in stagnation enthalpy, Cp(T t 2 ” T t j) 


AH* 

change in rothalpy, Cp (T* - T*) 


L 

Prandtl mixing length 


M 

Mach number 


P 

pressure 


r 

radius 


R 

radius of curvature, gas constant 


s 

entropy 


T 

temperature 


u 

relative velocity 


U 

cor, tangential velocity 


V 

absolute velocity 


W 

relative velocity 


y 

distance to nearest wall 


6 

boundary or shear layer thickness 


T1 

efficiency 



effective, laminar and turbulent viscosities 


V 

kinematic viscosity 


p 

density 


pu 

convective mass flow 


T 

shear stress 


CO 

angular velocity of the impeller 


Subscrij 

?ts 


abs 

absolute 


rel 

relative 


0 

tangential 
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Plane Subscripts 
1 flow inlet 

1 impeller inlet 

2 impeller exit 

3 vaneless diffuser at r^/^ ■ 1.0667 

State Subscripts 

0 reference stagnation, freestream at inlet 

t absolute stagnation 

tr relative stagnation 

u reached isentropically from i 

Efficiency Subscripts 

p polytropic 

s-s static-to-static 

t-t total-to-total 

w wheel 

Superscripts 

* rotary stagnation 

' reached isentropically 
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INTRODUCTION 


NASA Low-Speed Centrifugal Compressor Research Program [1] 

The new low-speed centrifugal compressor facility at NASA Lewis 
Research Center has been built to enable a more thorough understanding 
of subsonic flow in the complex geometric channels of centrifugal com- 
pressors to be obtained. NASA is planning an experimental program to 
obtain "benchmark" experimental data to verify three-dimensional viscous 
flow codes and to provide data with which to develop more sophisticated 
models of the various physical phenomena occuring in a centrifugal 
compressor. The complex phenomena to be considered include secondary 
flows generated by inlet vorticity, secondary flows generated when 
boundary layers are subjected to transverse forces, tip clearance flows, 
and turbulence suppression and enhancement due to strong normal pressure 
gradients occuring across the boundary layers in the impeller passage. 

NASA Low-Speed Centrifugal Compressor [1] 

The configuration of the compressor is shown in Fig. 1. The impel- 
ler has an exit diameter of 1.524 m (60 inches) and an inlet hub/tip 
ratio of 0.5. There are 20 full impeller blades with an exit backsweep 
of 55°. The design speed of 1920 RPM gives a rotor exit tip speed of 
153 m/s (503 ft/s). The design mass flow rate is 30 kg/s (66 lbm/s) at 
the standard inlet conditions of 101,325 N/m^ and 288.15 K. The rotor 
tip Reynolds number C^^^o) is ^ x Initially the compressor 
will be run with a vaneless diffuser. 

The nominal setting of the rotor tip clearance is a uniform 5.7 mm 
(0.225 inch). This tip clearance is 2.6 percent of the hub-shroud blade 
height at the impeller inlet and 4.0 percent of the blade height at the 
exit. 

The design philosophy adopted for the impeller was to obtain smooth 
velocity distributions while maintaining minimum blade lean. The pre- 
liminary design code used the Katsanis and McNally meridional flow 
program [2] for flow analysis. Velocity distributions were then calcu- 
lated using a three-dimensional Euler code as described by Denton in 
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reference 3. Figure 2 shows the velocity triangles and blade loading 
resulting from this preliminary design procedure. The rotor inlet tip 
Mach number is 0.31; the blade loading is quite uniform throughout the 
impeller; and the predicted relative Mach number at the impeller exit is 
apparently non-uniform, varying from 0.2 at the tip to 0.25 at the 
hub. The rotor exit absolute Mach number is 0.29. 

The NASA predictions of rotor total pressure ratio and (total-to- 

total) efficiency are then 1.173 and 0.934, respectively. The work 
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input parameter AH/U is anticipated to be 0.607. 

Contribution of Present Work 

The objective of this study was to perform a three-dimensional 
turbulent flow calculation for the NASA Compressor at its design operat- 
ing conditions. Similar calculations were recently performed by the 
authors for the earlier NACA 48-Inch Radial-Inlet Centrifugal Impeller 
[4,5]. Those calculations reproduced the measured overall performance 
and most of the features of the loss distributions in the NACA flow 
study. They gave further insight into the complex 3-D flow with tip 
leakage. It is therefore hoped that the present 3-D flow calculation 
for the NASA rotor will provide a useful prediction to aid in planning 
experimental work and flow measurements. 

For example, the preliminary 2-D boundary layer calculations of 
Wood et al. [1] suggested boundary layer thicknesses of approximately 
one inch (2.54 cm) near the NASA rotor exit. The present calculation 
shows regions of thicker boundary layer flow and loss accumulation where 
measurements may be easier, and also regions of thin boundary layer flow 
where measurements may be more difficult. In addition, the calculation 
includes a stationary shroud wall and a detailed representation of flow 
through the tip gap. This allows the shroud boundary layer development, 
including the infuence of tip leakage flow mixing, to be predicted; and 
the calculation shows regions of reverse flow in the meridional direc- 
tion. The development of the shroud boundary layer appears to strongly 
influence both the impeller and the diffuser performance. 
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The calculation results suggest flow mechanisms which are critical 
to the flow development and performance of this compressor. They dem- 
onstrate another level of preliminary flow calculation and prediction, 
for comparison with the results of the NASA preliminary design calcula- 
tions made for this compressor. The results also allow specific sugges- 
tion for changes in the impeller and diffuser geometries to improve the 
performance of the compressor to be made. 

COMPUTATIONAL METHOD 

The calculation was performed with the Moore Elliptic Flow Program, 
a 3-D pressure-correction solution procedure for discretized forms of 
the Navier-Stokes equations. The conservation equations for mass, mo- 
mentum and energy used for steady compressible flow are shown in Table 
1. The discretization of the conservation equations is described in 
reference 6. The equations are solved using the solution procedure 
described in reference 7. 

A Prandtl mixing-length turbulence model, summarized in Table 2, 
was used. For this calculation, a free-stream mixing length of 0.000387 
m was chosen, corresponding to an estimate of 0.5 percent free-stream 
turbulence at the rotor inlet. 

Use of Support Calculations 

Two preliminary 2-D flow calculations were performed: 

1. A meridional flow calculation was used to obtain inlet conditions 
for the 3-D flow calculation. 

2. A tip gap flow calculation was used to test the computational mesh 
chosen around the blade tip. 

COMPRESSOR GEOMETRY AND COMPUTATIONAL GRIDS 
Meridional Flow Calculation 

The inlet bellmouth, the inlet duct and the shroud contour of the 
impeller are shown in Fig. 3. 2-D meridional flow through the impeller 
was modelled using a 40 x 2 x 22 grid, shown in Fig. 4; this included 
the rotating spinner and hub of the impeller. 
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Table 1 . Steady Flow Conservation Equations 


Mass 

V*pu = 0 

Momentum 

pxi»7ii - (V«y7)ju « 7«yVu^ - 7p - 2p(jjxji - p^x(ojxr) 
Equation of State 

p * pRT 

Rothalpy (energy) 

pu'VH* - (V*yV)H* » 0 

H* = c T + l(u*u) - l(iuxr • ojxr) 

p ~2 T 

Table 2 . Prandtl Mixing Length Viscosity Model 


y = pL 2,, du" 
dy 

L is the smaller of 

0.08 times width of shear or boundary layer 
0.41 times distance to nearest wall 

Van Driest Correction 

L - 0.41 H y" (l.-expt-'yi/pT / 26 y i ] ) 

Near-Wall Correction 


3-D Flow Calculation 


The 53 x 21 x 22 grid for the 3-D flow calculation is shown in 
Figs. 5, 6, and 7. Figure 4 shows the inlet for the 3-D flow calcula- 
tion 0.188 m (7.4 inches) upstream of the impeller inlet, part way along 
the spinner; this is the start of the grids in Fig. 5. The calculation 
was performed through the vaneless diffuser to a radius ratio r/r 2 “ 
1.5. Figure 5a shows the basic constant-height vaneless diffuser; this 
was redesigned during the calculation to have a constant area starting 
at r/r 2 ■ 1.08, as seen in Fig. 5b. 

Blade-to-blade sections of the 3-D grid are shown in Fig. 6. The 
55 degree backsweep is evident in the shroud and hub sections in Figs. 
6a and 6b, and leading and trailing edge detail is seen in Figs. 6c and 
6d. Cross sections of the grid at the impeller exit and near the impel- 
ler inlet, shown in Fig. 7, illustrate the grid refinement used around 
the blade tip and the walls of the impeller passage. 

RESULTS OF MERIDIONAL FLOW CALCULATION 

In the 2-D meridional flow calculation, turbulent boundary layers 
were started on the shroud wall 1.04 m (41 inches) upstream of the 
location of the impeller leading edge and at the stagnation point on the 
spinner. The flow then developed as shown in Fig. 8. By the inlet for 
the 3-D flow calculation, the shroud wall boundary layer thickness was 
about 10 percent of the local passage height while the boundary layer on 
the rotating spinner was about 3 percent thick, as shown in Fig. 8c. 
The velocity vectors in Fig. 8a show separation of the shroud boundary 
layer near the impeller exit as the static pressure on the shroud wall 
recovers from the minimum at about 30 percent of meridional distance 
through the impeller, seen in Fig. 8b. At the exit of the flow domain, 
at r/r 2 » 1.5, the shroud wall boundary layer fills the diffuser but the 
negative mass flow is only about 0.7 percent of the total. 
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TEST OF GRID SELECTED FOR TIP GAP 


A two-dimensional tip gap/cavity geometry was set up to test the 
grid spacing within and around the blade tip. The ratio of the circum- 
ferential length to the height of the tip clearance gap varies from 
about one near the impeller leading edge to about two at the impeller 
exit. Figure 9 shows results for flow through a square tip gap of 
height h entering a cavity of overall height 10 h and width 9 h. A 
typical Reynolds number of 24,000, based on the mean velocity in the tip 
gap and the height h, was chosen. 

The velocity vectors in Figs. 9a and 9c also show the locations of 
the grid points. Overall, a 27 x 12 x 2 grid was used, with four grid 
points in the flow path between the blade tip and an inviscid shroud 
wall. The near-wall points are located 0.1 h from the blade tip; and 
with such a coarse grid there is no evidence of the expected flow sepa- 
ration at the sharp tip gap entrance. Indeed, although the flow dis- 
tribution at the tip gap exit in Fig. 9c seems qualitatively reasonable, 
it probably represents a tip leakage flow about 10 percent too large, 
due to an effective rounding of the inlet corner [8]. 

After leaving the tip gap, the flow mixes and spreads as it moves 
across the endwall, see Fig. 9b. It then splits as it encounters the 
cavity exit, with some of the flow circulating around the cavity. These 
same general features are probably superposed on the throughflow in the 
impeller. 

RESULTS OF 3-D FLOW CALCULATION 
Meridional Flow 

Straight Walled Vaneless Diffuser 

The 3-D flow calculation was started with the straight-walled 
vaneless diffuser. Figure 10a, however, shows that this resulted in 
strong backflow, in excess of 10 percent of the throughflow, at the 
diffuser exit. This violated the exit boundary condition. So the 
diffuser was redesigned to have a constant area from a radius ratio r/^ 
■ 1.08, i.e. from the leading edge of the proposed vaned diffuser, see 
Fig. 1. The area change was made by moving the hub wall, so that 
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changes in the development of the shroud wall boundary layer could be 
clearly seen. 

Modified Vaneless Diffuser 

Figure 10b shows the meridional flow development for the final 
converged solution with the modified vaneless diffuser. Now there is no 
backflow at the diffuser exit, although the shroud boundary layer is 
close to separation all through the diffuser. All the subsequent flow 
results presented here are for this modified diffuser. 

Flow near the pressure side of the impeller passage is shown in 
Fig. lObl. Velocity vectors for flow entering the tip gap are seen all 
along the shroud. Two other secondary flows also become evident in the 
second half of the passage; a secondary flow from the hub to the shroud 
due to curvature of the passage from the axial to the radial direction, 
and an opposing secondary flow due to the tip leakage flow impinging on 
the pressure surface. At mid-passage, Fig. 10b2 shows the shroud bound- 
ary layer thickening and separating in the impeller* 

Near the suction side of the impeller, centrifugal and curvature 
effects combine to give strong secondary flows towards the shroud. This 
is especially true in the second half of the passage, where flow angles 
as high as 30 degrees to the meridional direction are seen in Fig. 10b3. 

Thermodynamic Performance 

The thermodynamic performance of the compressor has been evaluated 
using mass-averaged properties from the 3-D viscous flow calculation. 
The procedure and nomenclature used follow closely those presented in 
reference 9. Table 3 presents the thermodynamic state points for the 
compressor at 1920 RPM and 30 kg/s. 

Four cross-sectional planes are used for these results. Plane 1 is 
the inlet to the flow domain, 0.188 m upstream of Plane i at the leading 
edge of the impeller blades. Plane 2 is at the impeller tip and Plane 3 
is a plane to be used in preliminary measurements in the vaneless dif- 
fuser at a radius ratio m 1.0667 (r = 32 inches). The states shown 
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are static states (1, i, 2, 3), total states (tl, ti, t2, t3), relative 
total states (trl, tri, tr2, tr3), and rotary stagnation states (*1, *i, 
*2, *3). Also shown is state u, reached by isentropic centrifugal 
pressure rise from state i, which is used to describe the impeller 
performance as a diffuser. The remaining states (t2* , tu’ , 3', t3’) are 
used in calculating the isentropic Impeller and compressor efficiencies 
in Table 4. The zero level for entropy is taken at the reference state 
( 0 ). 

The total pressure ratio at the rotor exit is 1.136 compared with 
the NASA prediction of 1.173 [1]. The static pressure ratio at the 
rotor exit is 1.087. 


13 


Table 3. Thermodynamic state points for the compressor at 1920 RPM and 
30 kg/s (66 lbm/s). 


State 

P/P 0 (1) 

T(K) 

s(J/kgK) 

0 

1.0000 

288.15 

0. 

1 

(0.9857) 

287.00 

0.1 

tl 

(0.9996) 

288.15 

0.1 

*1 

(0.9996) 

288.15 

0.1 

trl 

1.0263 

290.33 

0.1 

i 

(0.9786) 

286.40 

0.1 

ti 

(1.0004) 

288.20 

0.1 

*i 

(0.9998) 

288.15 

0.1 

u 

(1.0936) 

295.64 

0.1 

tri 

1.0296 

290.59 

0.1 

2 

(1.0867) 

295.95 

3.0 

t2 

(1.1362) 

299.74 (2) 

3.0 

*2 

(0.9911) 

288.26 

3.0 

tr2 

1.1388 

299.94 

3.0 

t2' 

(1.1362) 

298.87 

0.1 

2u' 

(1.0867) 

295.10 

0.1 

3 

(1.0931) 

296.62 

3.6 

t3 

(1.1348) 

299.82 (2) 

3.6 

*3 

(0.9893) 

288.28 

3.6 

tr3 

1.1584 

301.58 

3.6 

3' 

(1.0931) 

295.61 

0.1 

t3 ' 

(1.1348) 

298.79 

0.1 


^^Note: the slight inconsistency of these values arises because the 

pressures were independently mass averaged. 

^^ T t 3 should equal T^ for this adiabatic diffuser flow. 
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Table 4 


Efficiencies and work parameters for this compressor at 1920 
RPM and 30 kg/s (66 lbm/s). 


^s-slS 


n t-tl3 



- 0.895 


- 0.912 


pl3 


T T 

3 1 3 

in - — / in 

L l 1 


0.897 


w 


T - T 
2u» u 

T. - T 
2 u 


-1.74 


t-ti2 


T t2* T ti 
T t2 ' T ti 


- 0.925 


AH 

u 2 


c (I , - T ) 

p t2 tl 

2 2 
id r 2 


0.4960 


AH* 

u 2 


c (T* - T*) 
P . 1 L. 


2 2 
id r 2 


0.0047 


AH - AH* 




- 0.4913 


- 0.4914 
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Impeller and Compressor Efficiencies 

The total-to-total , static-to-static, and polytropic efficiencies 
for the compressor between states 1 and 3 are given in Table 4 as 91.2, 
89.5, and 89.7 percent, respectively. These are all significantly below 
the NASA prediction of 93.4 percent [1]. The total-to-total efficiency 
of the impeller between states i and 2 agrees more closely and is 92.5 
percent. These low values can be explained by the very low wheel effic- 
iency, n , which gives a measure of the rotor performance as a diffuser, 
reducing the relative kinetic energy. Table 4 shows this to be actually 
negative, at -174 percent (Note the small differences between large 
numbers here and the large errors possible in this figure. But it does 
seem to be negative!). 

Work Input 

2 

The work input parameter AH/U is seen in Table 4 to be 0.496. 

Most of this work is seen as a change in moment of momentum across the 

rotor, but some moment of momentum is dissipated at the shroud wall, 

between the inlet and exit of the impeller, and this appears as an 

increase in rothalpy (AH*) through the impeller. Thus the sum 

2 2 

of V^./air- - 0.4914 and AH*/U - 0.0047 should equal AH/U ■ 0.4960. 

The low value of work input is partly due to the high backsweep (55 
degrees) of the impeller blades. It is partly due also to blockage 
caused by the thick shroud boundary layer at the impeller exit, seen in 
Fig. 10b2. This causes high relative velocities in the hub half of the 
impeller passage, and these in turn cause reduced values of absolute 
velocity. The contours of dimensionless absolute tangential veloci- 
ty, Vq 2 /U 2 > seen in Fig. 11 in the exit plane of the impeller, therefore 
show a significant region with values less than 0.5 below passage mid- 
height. 

Static Temperature Changes in the Impeller 

The static temperature rise in the impeller may be expressed as 
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t 2 - t L 


2 2 2 2 
W - W U - U 

W i W 2 . U 2 U i 

T 


2c 


2c 


+ T* - T* 
2 L ± 


P P 

= (4.19 - 3.99) + (11.68 - 2.44) + (288.26 - 288.15) 


0.20 


9.24 


0.11 


m 9.55 K 


The centrifugal effects clearly dominate in this rotor and there is a 
surprisingly small contribution, only 2 percent, from the reduction of 
the relative kinetic energy. This is also seen in Fig. 12 where the 
development of the mass-averaged values of relative kinetic energy and 
kinetic energy associated with the local wheel speed may be followed. 
The relative kinetic energy falls linearly in the first half of the 
impeller passage but then it increases, due in part to excessive shroud 
boundary layer blockage and excessive backsweep. 


Total Temperature Changes in the Impeller 


The total temperature rise in the impeller may be written 


r t2 - T ti 


2 2 
W - W 
_i _2 

2c 


2 2 
U 2 ~ U i 
2c 


T* - T* 


2 2 
V - V 
V 2 V i 

2c 


P P P 

= (4.19-3.99) + (11.68-2.44) + (288.26-288.15) + (3.79-1.80) 


0.20 


9.24 


0.11 


1.99 


= 11.54 K 

Again the impeller performance is seen to be dominated by centrifugal 
effects. The absolute kinetic energy at the exit of the impeller is 
small, and there appears to be a limited potential for contributions 
from the diffuser downstream of the impeller. Said another way, the 
reaction of the impeller is large: in fact, 
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Reaction - ^ ^ - 0.828. 

3-D Pressure Distribution 

The three-dimensional distribution of reduced static pressure 
throughout the compressor flow domain is shown in Figs. 13, 14 and 15 as 
plots and contours of isentropic relative Mach number. 

Impeller Blade Loading 

The calculated blade loading for hub and tip sections of the impel- 
ler is shown in Fig. 13. Also shown, repeated from Fig. 2, are the NASA 
predictions made using a 3-D Euler code [3]. The present calculation 
confirms the uniform blade loading predicted by the Euler code, and the 
loading distribution at the hub agrees reasonably well. But there is a 
significant difference at the blade tip section. The 3-D flow calcula- 
tion with the pressure-correction method shows diffusion only in the 
first half of the impeller passage, followed by a slight reacceleration 
in the second half. As discussed above, this is partly due to shroud 
boundary layer blockage in the present calculation. Unlike the calcula- 
tion with the Euler code, the present results show a quite uniform rotor 
exit static pressure distribution. 

Blade Surfaces 

Figure 14 shows the surface Mach numbers on the pressure and suc- 
tion sides of the impeller blades, together with results on grid sur- 
faces extrapolated upstream and downstream. Figure 14a shows the mini- 
mum Mach number in the hub-pressure-side corner region to be about 0.17, 
and Figure 14b shows a maximum surface Mach number of 0.35 on the suc- 
tion surface in the inducer. Diffusion in the impeller is seen to be 
limited mostly to the shroud half of the passage and there it is com- 
plete by mid-way along the impeller. The static pressure is quite 
uniform from hub to shroud near the impeller exit and in the diffuser 
downstream. 
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Hub Wall and Blade Tip Region 

Isentropic Mach numbers on blade-to-blade planes are shown in Fig. 
15. At the hub, the Mach number varies from about 0.2 in the inducer to 
about 0.27 at the impeller exit, as seen in Fig. 15a. Near the shroud 
wall two sections are shown, one at about 80 percent of blade height in 
Fig. 15b and one at the height of the blade tip in Fig. 15c. The dif- 
ferences between the two planes are quite small in this calculation, of 
the order of 0.025 in Mach number in the inducer. Overall the Mach 
numbers change very little from about 0.31 in the inducer to about 0.28 
at the impeller exit. 

Shroud Static Pressure Distribution 

At the shroud wall, the static pressure distribution, normalized 
with the inlet total pressure, p Q , is as shown in Fig. 16. The pressure 
ratio varies from a minimum of .about 0.96 on the suction side in the 
inducer to an average value of 1.087 at the impeller exit. This steady 
pressure rise is mostly due to centrifugal effects. 

3-D Flow Development 
Blade-t o-Blade Planes 

Throughflow velocity vectors projected onto blade-to-blade planes 
are presented in Fig. 17. A plane near the hub, at a fraction C - 0.003 
of passage height, is shown in Fig. 17a, while Fig. 17b shows the re- 
sults at 80 percent of passage height. The flow appears quite uniform 
at the hub with little evidence of variations due to blade loading or 
impeller diffusion. At the 80 percent height, there appears to be some 
diffusion (slowing of the relative velocity) and some evidence of blade 
loading in the latter half of the passage, but these effects may also be 
due to a thickening of the shroud boundary layer. On both of these 
planes secondary flows appear small. 

Much stronger secondary flows are found at the mid-height of the 
blade tip gap, as seen in Fig. 17c. The convection of tip leakage flow 
across the shroud wall towards the pressure side can be seen in the 
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inducer; the leakage flow appears to reach the next tip gap about 40 
percent along the blade; and after that, there is backflow in the merid- 
ional direction in the shroud wall boundary layer. The tip leakage flow 
crosses the blade tip obliquely in the inducer and there is a strong 
component in the throughflow direction in the tip gap. In the latter 
half of the impeller, however, the tip leakage becomes more normal to 
the blade, contributing significantly to the backflow. 

Cross-Sectional Planes 

The non-orthogonal representation of secondary flow velocity vec- 
tors, introduced in reference 5, is used in Fig. 18 to aid in visualiz- 
ing the secondary flows. The results are shown on four planes, at 10, 
40, 74 and 100 percent of meridional distance, A, through the impel- 
ler. Contours of primary flow velocity, normalized by the impeller tip 
speed U 2 ■ 153 m/s, are also shown in an attempt to give a complete 
picture of the three-dimensional flow development. 

The strongest secondary flows are clearly associated with the tip 
leakage flow and its subsequent convection and mixing in the passage. 
As seen in Fig. 17c, the tip leakage fluid convects across the shroud 
wall; in Fig. 18a at 10 percent meridional distance it is just emerging 
from the tip gap; in Fig. 18b, at A - 0.40, it has just reached the 
pressure side and is flowing into the next tip gap. Some mixing and 
recirculation of the tip leakage fluid is evident in Fig. 18b at the 
edge of the shroud boundary layer at about 80 percent of passage height. 

In the second half of the passage, the tip leakage fluid impinges 
on the pressure surface of the next blade, causing secondary flow 
towards the hub. This was also seen in the meridional view of the 
velocity vectors near the pressure surface in Fig. lObl. By A ■ 0.74, 
the tip leakage flow and recirculation are well established features of 
the whole flow, as seen in Fig. 18c. Their influence remains at the 
impeller exit, seen in Fig. 18d; and they explain the distributions 
of V 00 /U 0 seen in the shroud half of the impeller passage in Fig. 11. 
The contours of 0.6 and 0.7 are caused by recirculating tip leakage 
fluid. 
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Of the other secondary flows evident in Fig. 18, the migration of 
the suction surface boundary layer towards the shroud wall is the most 
extensive. This was noticed particularly in Fig. 10b3. 

Detail of Flow through Tip Gap 

The calculated secondary velocity vectors and primary velocity 
contours for flow over the blade tip are shown in Fig. 19. Figure 19a 
shows the tip gap flow just starting in the inducer, while Fig. 19b 
shows a more developed leakage flow pattern with repeat leakage and 
entrainment of suction surface boundary layer flow. 

Cumulative Flow through Tip Gap 

The cumulative fraction of the passage mass flow that passed 
through the tip gap builds up quite linearly, as seen in Fig. 20. This 
is due to the uniform height of the tip gap and the uniform blade load- 
ing. The overall leakage is 18 percent, although most of the tip leak- 
age fluid appears to have passed over more than one blade tip. 

Loss Distributions 

Distributions of entropy on cross-sectional planes through the 
impeller and the vaneless diffuser are presented in Fig. 21. The con- 
tour interval is 1.0 J/kgK. 

In the inducer the flow is modelled as having turbulent boundary 

layers and these are of relatively uniform thickness. The suction 
surface boundary layer thickens by mid-passage and then convects towards 
the shroud, as seen in Figs. 21a, b and c. 

The accumulation of loss in the shroud boundary layer is the main 
feature of the impeller loss development. Initially the shroud boundary 
layer is thin (see Fig. 21a), but then it receives an input of tip 
leakage losses. The losses spread rapidly away from the shroud wall in 
the second half of the impeller passage (see Fig 21c and 21d). By the 
impeller exit, they cover most of the shroud/pressure-side quadrant of 
the passage. The rest of the flow is mostly low loss, potential flow, 
apart from a re-thickened suction side boundary layer. 
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The mixing downstream of the impeller is not very rapid. A compar- 
ison of Figs. 2 Id and 2 If shows some thickening of the hub and shroud 
boundary layers and a spreading of the blade wake. But the high loss 
region from the shroud/pressure-side quadrant is relatively unchanged by 
r/r£ “ 1.08. 

DISCUSSION OF RESULTS 
Tip Leakage Flow Mixing 

The NASA low-speed centrifugal compressor in its initial configura- 
tion has a tip clearance gap in the inducer much larger than the two 
impellers discussed in reference 9. The NASA rotor has a clearance gap 
equal to 2.6 percent of the inlet blade height, whereas in the DFVLR and 
Rolls-Royce impellers the gaps were 0.5 and 1.0 percent, respectively. 
The impeller has been set up for large tip leakage flows and that is 
what is found in the present calculation. The cumulative flow predicted 
to pass through the tip gap is 18 percent compared with 12 percent for 
the Rolls-Royce impeller. 

In passing over the blade tip, the fluid experiences increases in 
loss and turbulence kinetic energy. And, in the second half of the 
impeller passage, some of the leakage fluid passes through more than one 
tip gap. So there is enhanced turbulence in the boundary layer near the 
shroud wall. Then, too, some of this leakage fluid recirculates back 
towards the suction side after impinging on the pressure surface. So 
the outer layers of the shroud boundary layer may also have high turbu- 
lence. 

Offsetting this, for example, could be the curvature of the shroud 
wall which could interact with the boundary and mixing layers to sup- 
press turbulence. An estimate of the ratios, 6/R, of the shroud shear 
layer thicknesses to the radius of curvature of the streamlines suggests 
that these could be of the order of 0.1. Adams and Johnston [10] sug- 
gest that the Prandtl mixing length, L, could then be 0.0256 instead 

of 0.086 as modelled in the present calculations. The turbulent viscos- 

2 

ity which depends on L would then be a factor of ten smaller than 
predicted here. With the dominant influence of tip leakage flow and 
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mixing on the 3-D impeller flow development found here, it would seem 
useful to study the effect of such mechanisms on the turbulent mixing. 

Jet-Wake Mixing 

Jet-wake mixing in the vaneless diffuser at the exit of the NASA 
impeller appears to involve the jet in the hub half of the passage 
slowly mixing with high loss fluid in the shroud half. Superposed on 
this is a tangential spreading of the blade wakes which also appears to 
be quite slow in the present calculation. 

Backflow Along Impeller Shroud Wall 

Reverse flow in the meridional direction along the shroud wall is 
found from 33 to 114 percent of meridional distance in the impeller. It 
is seen in Fig. 17c as vectors with a negative meridional component, and 
it is seen in Figs. 18b, 18c, 18d and 19b bounded by contours (thick 
solid lines) of zero primary velocity. The picture in Fig. 18c at 74 
percent meridional distance is quite close to the section at 77 percent 
where the maximum backflow is about 2.5 percent of the throughflow. 
Here the thickness of the reverse flow region reaches a maximum of about 
1 cm. At this same impeller location, Fig. 19b shows that all the flow 
in the tip clearance gap has a component in the reverse direction. 

Backflow was also seen in the calculations by Moore and Moore [5] 
for the NACA 48-Inch Wheel. And its possible significance as a precurs- 
or of impeller stall has been discussed recently by Chen et al. [11]. 
Apparently such reverse flows may have a detailed vortex filament struc- 
ture. 

Comparison with NASA Flow and Performance Predictions 

Table 5 compares the performance predictions which have been made 
in this report with those made in the NASA design study [1]. Since the 
NASA study was for the compressor with its wedge-type vaned diffuser, 
comparisons are made for the rotor performance only. 
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The present calculation predicts a poorer performance for the 
current impeller configuration. The reasons for this different predic- 
tion are probably associated with the different impeller exit static 
pressure distributions and the different models of shroud boundary layer 
development and tip leakage flow mixing. The meridional flow and Euler 
codes used by NASA both predict higher impeller diffusion and corres- 
pondingly better compressor performance. 

Table 5. Comparison of Rotor Performance Predictions at 1920 RPM and 30 
kg/s (66 lbm/s). 



NASA [1] 

Present 

Calculation 

ft2 

1.173 

1.136 

1 ° <‘> 

0.934 

0.925 

t-tIZ 

AH 

0.607 

0.496 

u 2 

^reli ,tip 

0.31 

0.31 

M rel2 

0.20 

0.260 

^abs2 

0.287 

0.253 

Reaction 

0.763 

0.828 


(^Plane i is at the impeller inlet; Plane 2 is at the impeller exit. 
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CONCLUSIONS 


A three-dimensional turbulent flow calculation has been made for 
the NASA Low-Speed Centrifugal Compressor Impeller. The calculation was 
performed with the Moore Elliptic Flow Program, a 3-D pressure-correc- 
tion solution procedure for discretized forms of the Navier-Stokes 
equations. A Prandtl mixing-length turbulence model was used. 

The calculation was made at the compressor design conditions of 
1920 RPM and 30 kg/s, at 1 atm and 288.15 K. The tip clearance gap was 
a uniform 0.225 inches (5.7 mm). The constant-height vaneless diffuser 
was modified to have a constant area starting at a radius ratio r = 
1.08. 

The predicted performance is significantly worse than that pre- 

2 

dieted in the NASA study. The work done factor AH/U is 0.496 compared 
with 0.607. The total pressure ratio is also therefore lower, 1.136 
compared with 1.173. The total-to-total impeller efficiencies agree 
reasonably well, 92.5 percent in the present calculation and 93.4 per- 
cent predicted by NASA. 

The thermodynamic performance of the impeller is dominated by 
centrifugal effects. The efficiency of the impeller as a diffuser, 
reducing the relative kinetic energy, i.e. the wheel efficiency, is 
negative. The actual value calculated is -174 percent. Previous calcu- 
lations made by the authors for other centrifugal compressor impellers 
have given values of about +60 percent. But those impellers were run 
with significantly smaller tip clearances in the inducer, 0.5 and 1.0 
percent of blade height, compared with 2.6 percent, here. 

The tip clearance varied through the impeller from 2.6 percent to 
4.0 percent. Combined with the relatively uniform blade loading, the 
resultant tip leakage was 18 percent of the passage throughflow. The 
impeller flow development is therefore strongly influenced by the con- 
vection and mixing of tip leakage fluid. 

Tip leakage flow passes across the shroud wall; some of the flow 
impinges on the pressure surface of the next blade, while the rest 
passes through the next tip gap. A region of high losses develops in 
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the shroud/pressure-side quadrant of the impeller passage. This "wake" 
mixes relatively slowly with the "jet" in the hub half as the flow 
passes through the vaneless diffuser. 

The boundary layer on the shroud wall exhibits backflow in the 
meridional direction over the last two-thirds of the impeller length. 
The reverse flow is a maximum of 2.5 percent of the throughflow about 
three-quarters of the way along the impeller. Here, the maximum thick- 
ness of the backflow region was about 1 cm in the present calculation. 

The impeller with its initial uniform tip clearance setting will 
offer the opportunity to study a three-dimensional impeller flow domi- 
nated by tip leakage. Tip leakage mixing layers, similar to those 
predicted previously by the authors in a helicopter engine impeller, can 
be studied. In the NASA low-speed centrifugal compressor impeller these 
layers will be of the order of 1-2 inches thick. This should be an 
ample size for detailed turbulence measurements and flow definition. 

RECOMMENDATIONS FOR SECOND PHASE TESTING 

1. The geometry of the vaneless diffuser should be modified to prevent 
significant reverse flows on the shroud wall. In the present 3-D 
flow calculation, a constant area vaneless section, starting at a 
radius ratio r/r 2 ■ 1*08 (see Fig. 5b), accomplished this. 

2. Tip leakage flow mixing should be studied in the impeller. Partic- 
ular attention should be given to the shroud wall boundary layer and 
to loss accumulations in the shroud/pressure-side quadrant. Turbu- 
lence modification, for example due to curvature of the streamlines, 
may significantly affect turbulent mixing in these regions. 

RECOMMENDATIONS FOR THIRD PHASE TESTING 

1. The tip clearance gap in the inducer region of the impeller should 
be reduced to 1% of blade height or less. The reduction should be 
extended along the impeller until a wheel efficiency (see earlier 
discussion) of at least 0.5 is achieved. 
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Figure 1 . Configuration of the NASA low-speed centrifugal compressor. 
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Figure 3. Inlet bellmouth and shroud contours for low-speed centrifugal compressor 
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Figure 5 . Meridional view of grids for 3-D flow calculations. 
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Fig. 6d. Trailing edge detail 
at mid-blade-height 


Fig. 7a 



Fig. 7b. Near impeller 
inlet, at 4.1% 
meridional distance 
from leading edge. 

(r0,r) . 


exit (r9,z). 
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Fig. 8b. Static pressure contours, ^ ^inlet 
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Fig. 9cl. Static 
pressure contours 
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Fig. 9c2. Total 
pressure contours. 
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Fig. 9c3. Velocity 
vectors at the 
grid points. 
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Fig. 9c. Flow detail over simulated blade tip 





Fig. lObl. Near pressure side. 
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Fig. 10b. Constant area vaneless diffuser from r/r» = 1.08. 
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Figure 11 . Contours of dimensionless absolute tangential 
velocity, V 0 2 / U 2 » in the exit P lane of the 

impeller. Contour interval = 0.025. 

P, pressure side; S, suction side. 
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Figure 12. Contributions of centrifugal effects ( ) and 

relative kinetic energy ( ) to temperature 

changes and work in the impeller. 


48 




tOSWflN 
HDVW 3AI1V13M 


49 


Figure 13 ♦ Comparison of calculated blade loading with NASA predictions using 
a 3-D Euler code. 

A Calculated near blade tip, at 90% of passage height. 

• Calculated at hub. 



Figure 14 . Blade surface reduced static pressures 
represented as isentropic relative Mach 
numbers. Contour interval = 0.025 
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Figure 16. 


Shroud static pressure contours, p/p Q > over 
the impeller passage. 
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Fig. 17b. Near blade tip, C = 0.80. 
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Fig. 17c. Near mid-height of blade tip gap, C * 0.99. 
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Figure 19. Calculated secondary velocity vectors and primary 
velocity contours for flow through the tip 
clearance gap. 
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Fig. 21d. A « 1.0, impeller exit. 
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